clear;
clc;
A=[16,2,3,13;
    5,11,10,8;
    9,7,6,12;
    4,14,15,1];
Q=orth(A),
norm(Q'*Q-eye(3))